Grassmann tensor network states and its renormalization 
for strongly correlated fermionic and bosonic states 
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The projective construction (the slave-particle approach) has played an very important role in un- 
derstanding strongly correlated systems, such as the emergence of fermions, anyons, and gauge the- 
ory in quantum spin liquids and quantum Hall states. Recently, fermionic Projected Entangled Pair 
States (fPEPS) have been introduced to efficiently represent many-body fermionic states. In this pa- 
per, we show that the strongly correlated bosonic/fermionic states obtained both from the projective 
construction and the fPEPS approach can be represented systematically as Grassmann tensor prod- 
uct states. This construction can also be applied to all other tensor network states approaches. The 
Grassmann tensor product states allow us to encode many-body bosonic/fermionic states efficiently 
with a polynomial number of parameters. We also generalize the tensor-entanglement renormal- 
ization group (TERG) method for complex tensor networks to Grassmann tensor networks. This 
allows us to approximate the norm and average local operators of Grassmann tensor product states 
in polynomial time, and hence leads to a variational approach for describing strongly correlated 
bosonic/fermionic systems in higher dimensions. 
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I. INTRODUCTION 

Traditional condensed matter physics is based on two 
theories: symmetry breaking theory for phases and phase 
transitions, and Fermi liquid theory for metals. Within 
the Fermi liquid theory, one assumes that the ground 
state wave function for the electrons can be approx- 
imately described by a Slater determinant. In other 
words, one assumes that the many-electron ground state 
can be constructed by filling the single-particle energy 
levels. Such an energy-level filling picture becomes a 
foundation for traditional many-body physics. In this pa- 
per, we will call states obtained by filling single-particle 
energy levels as energy-level-filling (ELF) states. In ad- 
dition to Slater determinant states, fermion paired states 
are also ELF states. 

We may view the ELF construction of many-body 
states as an encoding method of a physically relevant sub- 
set of states. Although a random many-body state can 
only specified by an exponential amount of data, hence 
making it impossible to specify and calculate physical 
properties efficiently, physically relevant states seem to 
have a much simpler entanglement structure. A generic 
ELF state on a lattice can be written as 



'J'/) = exp 



' J ij C ] C i 



|o> 



(1) 



where Y]uj\ sums over the pairs of sites in the lattice. 
Here we consider a simple example of paired spinless 
fermions, and unpaired fermions can also be represented 
as limits of such states. Such a many-body fermionic 
state is encoded by polynomial amount information char- 
acterized by Uij. A crucial property is that it is very 
easy to calculate the norm and the averages of any local 
operators for a ELF state; here "easy" means that the 
computational cost scales as a polynomial in the num- 



ber of modes. This effective encoding and the ease of 
calculating physical quantities (such as energy) form the 
foundation of the standard mean- field theory for interact- 
ing electron systems. In such a ELF approach, we may 
view as variational parameters and minimize the av- 
erage energy by varying uy. The Uij that minimize the 
energy give us an approximated many-fermion ground 
state. From the form of Uij we can determine which 
symmetry is broken and obtain the phase diagram of the 
system. 

However, the ground states for some strongly corre- 
lated electron systems cannot be approximated by ELF 
states, i.e. they cannot be constructed by filling energy 
levels. One classic example is the filling fraction v = 1/3 
Laughlin state 1 



<F, 



IK' 



(2) 



Although the v = 1 state *i = fT* 



^3 



li<j( Z i "]) 

is a ELF state, its cubic power, ^> 3 ~ (^i) 3 j is very 
different from any ELF states. In order to obtain the 
low energy effective theories for systems that cannot be 
described by ELF picture such as spin liquids and non- 
Fermi-liquid metallic states, a projective construction 
(also known as slave-particle approach or parton con- 
struction; see appendix B for a brief introduction) was 
developed 2-12 . Those states may appear in high T c super- 
conductors and other strongly interacting systems, and 
are by now widely used. Just like the ELF states, the 
projective states can also be characterized by a polyno- 
mial amount of information. So the projective approach 
can also be viewed as an efficient encoding of many-body 
states. 

In the projective approach, we view the projective 
states as variational trial wave functions 9 ' 11 for obtain- 
ing approximate ground state. What makes the projec- 
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tive approach so attractive is that from the form of the 
projective state, we can usually obtain the low energy 
effective theory that describes low energy excitations 8,11 . 
From the low energy effective theory, we have learned 
that projective states can capture many qualitatively 
new phenomena that ELF states fail to describe, such 
as fractional charge, fractional statistics and topological 
orders 7 ' 8 ' 12 . 

However, it is much harder to calculate the norm and 
local expectation values in the projective approach. Al- 
though expectation values can efficiently be calculated 
using variational Monte Carlo methods in the case that 
electron or spin operators can be expressed as the prod- 
ucts of parton operators 13,14 , Monte Carlo fails when 
those operators are expressed as sums of the products 
of parton operators 8 ' 15 ' 16 . As such there is no general ef- 
ficient way to calculate the norm and the average of local 
operators. 

In a different development, concepts of quantum in- 
formation theory have allowed one to gain a better in- 
sight into the entanglement structure in ground states 
of generic Hamiltonians of strongly correlated quantum 
spins. It has been shown that ground states of local 
Hamiltonians obey so-called area laws, and that ground 
states can therefore be efficiently be represented by the 
class of so-called Matrix Product States in one dimen- 
sion or Projected Entangled Pair States (PEPS) in higher 
dimensions 17 ' 18 . This class of states share some resem- 
blances with the states obtained in the projective ap- 
proach, as they can be understood as projections of lo- 
cally maximally entangled pairs. Importantly, techniques 
have been developed that make it possible to efficiently 
calculate expectation values of local operators for this 
class of states, as this can be done by contracting net- 
works of tensor products. In this paper, we will show 
how to generalize such networks to include Grassmann 
tensors, and show how they can be contracted efficiently. 
This allows us to define a large and important subclass of 
the states obtained in the projective approach for which 
it is possible to calculate expectation values. It also al- 
lows us to generalize all tensor product state methods to 
the fcrmionic case; special subclasses include the recently 
introduced fcrmionic PEPS (fPEPS) 19 and the fermionic 
Multiscale Entanglement Renormalization Ansatz 20 . 

The Grassmann tensor product states allow us to ex- 
press the norms of the projective wave functions in terms 
of Grassmann valued tensor network. Similarly, the av- 
erage of any local operators for the projective wave func- 
tions can also be expressed in terms of Grassmann tensor 
networks with a few "impurity" tensors. 

Throughout our study, we will also find that the Grass- 
mann tensor product states can be more general than 
projective states. Using Grassman tensor networks, we 
can systematically construct more general strongly corre- 
lated states for both bosonic and fermionic systems. We 
only need polynomial amount of data to characterize the 
tensor network. Thus the Grassmann tensor network ap- 
proach gives us an effective encoding for both fermionic 




FIG. 1: A random lattice where the physical degrees of free- 
dom are localized on the vertices. 



and bosonic many-body states. 

Calculating the norm and the expectation value of 
local operators for a tensor product state can be ex- 
ponentially hard in general 21 . However, many possible 
polynomial approximation schemes, including the tensor- 
entanglement renormalization group (TERG) method, 
have been proposed in recent years 22-27 . In this paper, 
we will show how the TERG method 24 can also be ap- 
plied to Grassmann tensor network. In particular, this 
implies that the average energy and other physical quan- 
tities of a Grassmann tensor network state in two dimen- 
sions can be calculated efficiently by using tensor net- 
work renormalization approach. This Grassmann tensor 
network approach is the natural approach for express- 
ing fermions in tensor network methods and hence pro- 
vides a new starting point for studying strongly corre- 
lated bosonic/fermionic systems. 



II. TENSOR-NETWORK REPRESENTATION 
OF ELF STATES 

Let us first consider a simple example of using tensor 
networks to represent a ELF state on an arbitrary graph, 
where the fermions live on the vertices of the graph (see 
Fig. I): 



|*/) = exp 



As 



io)=n( i +^ c H)io> (3) 



-MjjctcJ, W.l.O.g. Uij 



v^j^i — — "tj "-j "-j i w.i.u.5. u,jj — —Uji. The many- 
body wave function \&/({mj}) is given by 

= <oi n^) m< n i 1 + wki) i°) ( 4 ) 



where n%i = 0, 1 indicating if the site-z is empty or occu- 
pied. Note that the fermion operators Cj in the product 
n,( c i) mi 1S ordered in the following way 



H^r- =(ci) mi (c2) m2 (c 3 ) m3 ... 



(5) 
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Motivated by fermionic path integral, we may repre- 
sent the above wave function in terms of Grassmann num- 
bers Oi and their derivatives d6i, which satisfy 



0i6j — — OjOij d9i ddj — — dOjddi, 

dOil = 0. (6) 

We find that the wave function can be rewritten as 



dOidj — Sij 



T} = Mi, T? = 1, Gn = 1 + UijOjOi, (7) 

where J "integrates out" all Grassmann numbers. 
Similarly, ^({m^}) can be expressed as 

*/(«)= /n^n( i+u w)> 

i ij 



Tl = d0 u T° = 1, 

Gij = l-</^ = l + </ i 7 , (8) 
Note that J\ i and have different orders: 

i i 

(9) 

Here we have used the following identity 

Thus the norm of the wave function is given by 

<*/i*/> = e An^rn^xn^n^) 

{m,} i ij « ij 

(11) 



=/n^n^ 

where 

Tj = 1 + MiMi, = (1 + </^)(l + UijOjOi). 

(12) 

We may view as a dimension- 1 tensor on the vertex 
i, and Gij as a dimension- 1 rank-2 tensor on the link ij. 
Then f\ i T™ ; FJ^ Gij can be viewed as the tensor trace 
on such a tensor network. Note that the tensors contain 
Grassmann numbers and TVs always appear in front of 
Gi 3 's. We see that the norm of a fcrmion wave func- 
tion can be expressed as the tensor trace of a Grassmann 
tensor network. 




FIG. 2: A strongly correlated system on a graph. The physi- 
cal degrees of freedom live on the vertices. The vertices (the 
filled dots) are labeled by i, j, etc . The open dots are la- 
beled by ij, etc . The links connecting open and filled dots 
are labeled by /, J, K, etc . 



III. GRASSMANN TENSOR-NETWORK 
REPRESENTATION OF GENERIC STRONGLY 
CORRELATED FERMIONIC AND BOSONIC 
STATES 



The above result can be generalized to strongly cor- 
related fermionic states as well as strongly correlated 
bosonic states. Let us assume the physical degrees of 
freedom are localized on the vertices of a graph (see Fig. 
1). The vertices are labeled by i, j, etc . The states 
on vertex i are label by mi = 1,2,... Let a sign func- 
tion s(mj) to have the following property. s(m,i) = 1 if 
the state mj is a bosonic state and s(rrii) = — 1 if the 
state rrii is a fermionic state. To construct a many-body 
wave function \I>({mj}), we introduce some fermion op- 
erators ip? on each vertex i with a denotes the fermion 
species (such as spin). We also use ij etc to label the open 
dots and /, J, etc to label the links between the open and 
filled dots in Fig. 2. Then we can construct ^({rm}) as 



*(W)=E(°in c '^« 

{ ai } % 



(ij) 



|0>, (13) 



where K, L, etc in C mi - aK a L ... label the links that con- 
nect to the vertex i, and /, J in tij- aiaj label the links 
that connect the vertex i, j with the open dot ij. All 
the link indices a,K,aL,ai,aj ■■■ — 1, ...,.D are the inner 
indices defined in the standard (bosonic) tesor product 
states(TPS). Here tij- ai bj only contains Vf 1 ^ an d ip"^ 
operators and C mi - aK a L ... only contains tp" operators. 
tij-,aiaj always contains even numbers of fermion opera- 
tors. C mi - aK a L ... contains an even numbers of fermion op- 
erators if s(ffii) = 1 and contains odd numbers of fermion 
operators if s(m,i) = — 1. 

Introducing the Grassmann numbers Of for each ver- 
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FIG. 3: A tensor network formed by two kinds of Grassmann 



tensors T- 



and Gi 



etc 



tices) are labeled by i which represent T t .J 



The filled dots (the ver- 
, and the open 

dots are labeled by ij which represent Gij- aiaj - The lines (la- 
beled by /, J, etc ) connecting the dots represent the indices 
ai of the tensors. For two tenors connected by a line, the 
values of the associated indices in the two tensor are set to be 
equal, and those values are summed over in the tensor trace. 
Such a summation is represented by m 0-^)- 



tex, we can express the many-body wave function as 



*({mj) 



(14) 



where T™ l KQi is obtained from C mi - aKaL ... by replac- 
ing if>f by dOf and Gij- aiaj is obtained from % ;ai oj by 
replacing ipf* by Of. So, \&({mj}) can be expressed as a 
tensor trace over a Grassmann tensor network (see Fig. 
3). 

Similarly f can be expressed as 



= E 

{ai} 



Yi T t?a K a L ...Y[ G 



y;a/Oj ; 



(15) 



r.., r(aL ... is obtained from C^..^^... by replac- 



, r . aiaj is obtained from t\.. aiaj 



where T i 

ing Vf f by d0f and G 
by replacing ?/>f by Of. Essentially, (14) is a fermionic 
generalization of the standard (bosonic) tensor prod- 
uct states(TPS). If there is no Grassmann number in 
T, G, it becomes the standard TPS (one can further put 
Gij- ai aj = fiai,aj if the tensor contraction of inner indices 
are made over a trivial metric). 

The norm of the wave function can be calculated in 




FIG. 4: the Grassmann number Of is associated with the link 
/ that connects to vertex i. The tensor Gij contains only 
Grassmann numbers Of and 0". 



the same way: 

<*i*> = e / (n^u..n e ^^) x 

{mi}, {a/a/} i ij 



/ n^joVoi'-^ijaKOi" II G ir,aiaj G ir,aiaj 
' ij 

J^Gij-^jaj ,ajaj (16) 



{m i },{a I a I } i 



{a,a I } i ij 



where 



rp \ ^ rpm rpra 

J- i;a,K<iK -aL^L ■■■ / , 1 i\aK&L--- i;aKCiL--- 



Grij^ajai^ajaj — Gij^ aiaj Gij^ aia j 



(17) 



Again, the norm is a tensor trace of a Grassmann tensor 
network (see Fig. 3 where T, G are replaced by T, G and 
each link is indexed by a pair aj a/). We may combine 
the pair of indices (a/a/) into one pi, and rewrite the 
above as: 



(18) 



{Pi} » 



IV. 



GRASSMAN TENSOR PRODUCT STATES 
A "BOND" FORM 



Although the Grassmann tensor network representa- 
tions for the norm of fermion wave function, (11) and 
(18), are simple and compact, however, it is not easy 
to implement the renormalization calculation for such a 
generic Grassmann tensor network. In this section, we 
will consider Grassmann tensor networks that have a spe- 
cial form which makes the renormalization calculation 
easier. Although the Grassmann tensor networks have 
a special form, they can still represent generic strongly 
correlated fermionic states described by the more generic 
form (18). 

In the previous Grassmann tensor network, the Grass- 
mann numbers on vertex i, Of' , are labeled by a* = 



■5 



1,2,... To construct the more special form of the Grass- 
mann tensor network, we separate the Grassmann num- 
bers on vertex i into several groups, one group for each 
link that connects the vertex i. So each group of Grass- 
mann numbers is actually associated with a link labeled 
by I. Thus it is more convenient to relabel the Grass- 
mann numbers on vertex i as Of , which correspond to 
the Grassmann numbers on vertex i and associated with 
the link I that connects to the vertex (see Fig. 4). 

With the new labeling scheme, we can specify the spe- 
cial form of the Grassmann tensor network by requiring 
Gij- aiaj to only contain the Grassmann numbers Of and 
8j (see Fig. 4). (Note that Gij ;aiaj still contain even 
numbers of Of and Ofj.) Under such scheme, the wave- 
function (14) can be represented as: 



{aj} i 



ly.azaj , 



(19) 



where 

i;OKOt... / j *-i;a K ai J .. I > 

{i a K K Wl L }- Ieta ' 

(20) 

Here I € i means the links that connect to vertex i. 
Note that, in the expression of Gij, we have assumed 
that the link / connects vertex i with the open dot ij, 
and the link J connects vertex j with the open dot ij (see 
Fig. 4). We also note that Iff — or 1 indicates the 
presence or the absence of the Grassmann number dOff . 
So we may interpret Iff as an "occupation number of 
fermions" .Here ^2n a K }{i" L } sums ovcr & U the possible 
"occupation" distributions {Iff = 0, l\ax = 1,2,...}, 
{lf L = 0, l\a L = 1,2,...},... Similarly as in (14), 
S/ei ~ represents a fermionic state and 

Yliei J2 ai tf 1 ~ cven represents a bosonic state. How- 
ever, J2 ai If +J2 a , Kj' snou ld always to be even. Such 
type of representations for fermionic wave functions was 
first introduced in Ref. 19. 

In the new form, the norm of wave function can be 
expressed in the same way as in (18) with: 



i\a K a K ,a L a L . 



G 



ij;aiai 



ajaj Gij- aia jGij- aia j . 



(21) 



Here T™ &k&l and Gij- Sl aj ar e defined as: 



E 



n m i; {^}{^}. 



lei a r 



E Wrf'IKor) 

{??'}{?/ 7 } a ' aj 

(22) 



By combining the pair of indices (ajaj) into one pi 
and reordering those d6f',d6f'{6f',6f'), the tensors T 
and G in (18) can be further expanded as 



T 



T {n a K K }{r 

i;PKPL- 



E 

{n a K K }{n a L L }.. I&a 

Ga, PIPJ = e $&if jJ} nw j r" u^T' 1 

{«?}{«/'} 

(23) 

where the group of Grassmann numbers {rff'\oti = 
1,2,...} is the combination of {6j\j3i = 1,2,...} and 



{0f 1/3/ = 1,2,...}. Notice T^ kK}{t 



:\PkPl- 



vanishes if the 

total numbers of Grassmann numbers drff is odd and 
,K J }{"" J } 

^ij-.PlPJ 

numbers T]f,T]j is odd: 



{ n K K }{ n L L }• 



PkPl- 



= 0, if EE< J = odd < ( 24 ) 



^-° J} = o, if 5>?'+$>s'= odd. 



To gain some intuitive understanding of the new form 
of the Grassmann tensor network, let us represent the 
free fermion wave function \& f (4) using the new form of 
Grassmann tensor network: 



i ij 

Tt = J2<iei, T? = l, c 



l + u 



ijVjVl, 



(25) 



lei 



where J2 Iei sums over all the links that connect to vertex 
i. The norm of such a free fermion state is given by 

<*/i*/>= / n^n^ 

i ij 

T, = l + E d^][E 



Kei 



Lei 



Gij = (1 + <//0j)(l + UijOjOj). (26) 

Again, in the expression of Gij, we have assumed that 
the link I connects the vertex i with the open dot ij, and 
the link J connects the vertex j with the open dot ij. 

We would like to stress that the above formulism is 
very general. It can be used to express free fermion 
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states, as well as strongly correlated fermionic/bosonic 
states obtained from the projective construction. It can 
even express strongly correlated fermionic/bosonic states 
beyond the projective construction. We also would like 
to stress that the graphs discussed in this paper do not 
have to be a real space lattice. 



V. STRONGLY CORRELATED STATES FROM 
THE PROJECTIVE CONSTRUCTION 

As we have mentioned, the Grassmann tensor net- 
work can represent very general strongly correlated 
fermionic/bosonic states. In this section, we will con- 
centrate on strongly correlated fermionic/bosonic states 
obtained from a projective construction and show that 
all of these states can be expressed as Grassmann tensor 
product states. The Grassmann tensor network for such 
projective states takes a particularly simple form. 

To construct the Grassmann tensor network for the 
projective states, let us first construct the Grassmann 
tensor network for a ELF state: 



| "J/) = exp 



(ij) 



|0> 



(27) 



where there can be several kinds of fermions labeled by 
a, /3 on each vertex. As discussed in section IV, the above 
ELF state can be expressed in terms of Grassmann tensor 
product state: 

*/({m ilQ }) = (0|niI(^) mi ' a l*/) 

i a 

i ij 



T 



{ mi , a } 



a lei 



G, 



exp 



E 

(ij) 



(28) 



where 

„th 



0, 1 describes the occupation of the 

■,{ro», a } 



a"" fermion on vertex i. 1 is obtained from 

Y[ a {'4>t,a) mi - a by replacing by J2iei dOf where J2i ei 
sums over all the links that connects to the vertex i. Ga 



is obtained from exp 



^ {ij) {u l3 ) a ^]^l a \ by replacing 

ipjpiPia by ^j^j where / and J label the links that con- 
nects to the open dot ij between the two vertices i and 
j. Note that (28) is a generalization of (25). 

To obtain a strongly correlated fermionic state from 
the above ELF state, for example, we may assume a = 
1, 2, 3 and consider the following many-body wave func- 
tion: 



*corr(W) = (0| n^r*!*/), Ci = I| ^ 



(29) 



where mi — 0, 1. We note that v I / corr({'Tii}) is a many- 
body wave function for a spinlcss fermion system. It is a 
strongly correlated projective state which is very different 
from any ELF state. Such a state can be expressed in 
terms of Grassmann tensor network 



*corr(W)= J UV'U^iJ 



(30) 



nzO' 



G t 



exp 



a=l lei 



(ij) 



Similarly, to obtain a strongly correlated hardcore 
bosonic state (such as spin liquid state) using projective 
construction, for example, we may assume a = 1, 2 and 
consider the following many-body wave function: 

* S pin(W) = (0| n^)" 1 *!*/), h = Mi,2, (31) 



where m; = 0, 1. v I'spin({'7ii}) is a many-body wave func- 
tion for a hardcore boson system. It can also be viewed 
as a wave function for a spin-1/2 system. Such a state 
can be expressed in terms of Grassmann tensor network 



*»pin({m i })= / 1 1 // \\G, r 

i ij 

2 

T r = (IlE d ^) m, < g ^ = cx p 



(32) 



a=i lei 



^ijU^ef 

(ij) 



Both strongly correlated states ^co^dmi}) and 
^spin({wi}) are parameterized by (uij) a p. We may view 
{uij)aj3 as variational parameters. After finding (uij) a /3 
that minimize the average energy, we obtain the approx- 
imated ground state. We can also obtain the low energy 
effective theory from the form of the ansatz (uij) a p. 

Finally we would like to point out that the ways to 
construct ELF states and projective states discussed in 
sections II and V are simple but not efficient. The Grass- 
mann tensor product states derived in that way will usu- 
ally contain long range connections, which is not neces- 
sary in special cases 19 . 



VI. FERMION COHERENT STATE 
REPRESENTATION 

In above sections, we have represented the Grass- 
mann tensor network wavefunctions under the Fock ba- 
sis. Although the Fock basis representation is simple 
and straightforward to derive, however, because of the 
anticommutating relations for different Grassmann num- 
bers, the wavefunctions depend on the ordering the lo- 
cal Grassmann tensors T™,* „ and are inconvenient for 
simulating physical quantities for fermion systems. 

In this section, we would like to introduce the fermion 
coherent state representation for Grassmann tensor prod- 
uct states. We show the fermion wavefunctions in this 
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basis are independent of the ordering of local Grassmann 
tensors. To see this explicitly, let us consider a simple 
spinless fermion tensor product state: 

(33) 

where = 0, 1 represents the fermion occupation num- 
bers. It is easy to check that we can derive the Grass- 
mann tensor product wavefunction (19) under the follow- 
ing Fock basis: 

n(cirio) = (ctr i (4r 2 (4r 3 -io> (34) 

i 

The over complete fermion coherent state basis is de- 
fined as: 

1^=11(1-^)10), (35) 

i 

with closure relation: 

fl[d^drH(l-v:rh)\vHv\ = 1- (36) 

It is easy to derive the wavefunction for Grassmann ten- 
sor product state (33) under such a basis: 



then the wavefunction can be represented as: 



{vm}, {a/} i ij 



Notice under the new definition, T™? „ always contain 
even number of Grassmann numbers hence the definition 
of the wavefunction (39) is independent of how we order 
those local Grassmann tensors TJ"* „ 

It is very convenient to use the above wavefunction 
to calculate local physical quantities. The only thing we 
need to take care is the over complete nature of the basis, 
hence a proper measure is needed when we calculate the 
inner product for a wavefunction. For example, the norm 
of the wavefunction (39) can be calculated as: 



= EZ J IT Ti; PKPL ... Y[ Gijwipj ! (40) 

{pi} i ij 



w=ee/ n^r^ ,...n G ™- (37) 

{mi} {ai} i ij 

If we redefine the Grassmann tensor TJ™* „ as: 

{lK K Wl L }- Ieta ' 

(38) 



which is exactly the same result as we derived be- 
fore (tensors T and G are defined in Eq. (21)). 

Other physical quantities like the average energy can 
also be easily calculated in a similar way. For any local 
operators containing and c, we only need to replace 

with rj* and c with 77, then integrate out the Grass- 
-mann number respect to a proper measure as we do for 
calculating the norm. For example, the nearest neighbor 
paring term cjct can be expressed as: 



J i' 

EV- / rpl rpl /=( _ _ rpQ rpQ (-1 

/ j I 1 j\a M a,N ■■ i\a, K a L ..^ij;aiaj ±i:a K a L .A j\a M aN ■■ i j'- a i a -> 



EZ y Ti; aK a K ,a L a L ..Tj ;aM a M ,a N a N ...Gij;a I a I ,ajaj Y\_ ^ i' ;a K i a K , ,a L ,a L , ... H G%' j'\apap ,a ji a ji , (41) 
{a pap} V i'j' 



where the impurity tensors T 1 and T" are defined as: and the indices i' denote other sites beside i and j. 

rpl (Til rpO 

i\a K a K ,a L a L ... i;a K a L ... i;a K a L ... 

T" - - — f 1 - - T° (42) 
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In the last line we omit the minus sign because the 
two impurity tensors T and T contain odd number of 
Grassmann numbers hence they anticommute with each 
other. 

In conclusion, the norm of a Grassmann tensor prod- 
uct state can be expressed as a tensor trace of uniform 
Grasmann tensor-net and other local physical quantities 
such as energy can be expressed as a tensor trace of Gras- 
mann tensor-nets with impurity tensors. Here we have 
already seen that the fermion coherent state representa- 
tion is extremely convenient for expressing the norm and 
physical quantities of Grassmann tensor product states as 
tensor traces of Grassmann tensor-nets. Essentially, such 
representations are the most natural representations for 
Grassmann tensor product state and provide us a deep 
insight into what is a fermion wavefunction. Actually, 
a fermion wavefunctions should be described Grassmann 
numbers rather than complex numbers, a detail discus- 
sions will be represented in our future publications. 



trace over a Grassmann tensor network is a very formal 
exercise. Calculating such a tensor trace directly is a ex- 
ponential hard problem. So it seems that we gain nothing 
from writing the norm in the form of tensor trace. 

However, if we only want to evaluate the norm and 
physical quantities approximately, then there is a poly- 
nomial way to do so in terms of the Grassmann tensor 
network. The basic idea is to perform a coarse graining 
transformation of the tensor network which can simplify 
the tensor network into one with only a few tensors. 22,24 

In this section, we will explain how to apply a coarse 
graining transformation to a Grassmann tensor network. 
We will discuss three basic moves on a honeycomb lat- 
tice, using rank-two, rank-three and rank-four tensors as 
examples. 

In the first move, we combine two rank-three tensors 
Ti, Tj and one rank- two tensor Gij into one rank-four 
tensor T (see Fig. 5). We have 



VII. THE COARSE GRAINING 
TRANSFORMATION OF GRASSMANN TENSOR 
NETWORK 



P1P2P3P4 



/ ^i;p5PiP2^j;p6P3P4^' ! ij;p5P6 (43) 

P5P6 Ji l 



Expressing the norm and physical quantities of a 
strongly correlated wave function in terms of a tensor 



J 



where J.j only integrate over ry" 5 and ?y" 6 . If we expand 
Ti, Tj, G^, and T u we get 



■ j;P6P3P4 



e 4fp 1 pf 1}{ "" } n( d % 5 ) <5 n( d ^ i ) nri n( d ^ 2 )" 2Q 



\p r K 6 H«3 3 }{«4 Q4 } 
2-^i i;p6P3P4 



*6 t\"3 



P1P2P3P4 — 

K 5 >K 6 } 



Gij;p 5 p 6 — G, 
K 5 H< 6 } 



u;P5P6 



n(do <6 n( d % a3 )" 3 " 3 n( d ^ 4 ) <4 

a e a 3 a* 

ai CK2 Ct3 a 4 

n« 5 )" ?5 n« 6 ) <6 w 



We find that 



'P1P2P3P4 — 2^1 ill 



•\P5PlP2 



^ij;P5P6 'j:,P6P3P4 



(45) 



p 5 P6 k^k 6 } 



This allows us to calculate T from T t , Tj, and Gij. 



method discussed in Ref. 22) 



The second basic move splits the rank-four tensor T 



into two rank-three tensors T 



Tj and a rank-two tensor 



6). We first rewrite Tp^plplp™ 2 



}{«3 3 H< 4 } 



dj (see Fig. 

as (say, using the singular-value-decomposition(SVD) 



% 



.{n 2 * 2 }{™ 2 " 2 }« 3 }K 4 } 



i;qpiP2 j;qp3P4 



(46) 
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The above decomposition can be rewritten as 



T 



T-{«2 2 }{™2 2 }{' 
•PW2P3P4 



! }{< 4 } 



E 

P5P6 



MP5P1P2 y ij;p: 



T,P6P3P4. 



(47) 



FIG. 5: We note that the Grassmann numbers 6"' and d6"' 
are associated with the links the connect the G tensors and 
T tensors. Here we label the links with J = 1, 2, 3, 4, 5, 
and 6. The tensor Ti- P5P1P2 contains Grassmann numbers 
{d??5 5 , dr/^ 1 , d??2 2 }■ The tensor Tj- P6P3P4 contains Grass- 
mann numbers {dr/Q 6 , d^g 3 , dr?J 4 }. The tensor T contains 
Grassmann numbers {drj" 1 , drj^ 2 , dr)^ 3 , d??" 4 }. The tensor 
dj contains Grassmann numbers {% 5 ,r/Q 6 }. 




where 



and 



->{n 5 }{n e } 

^ij;pnP6 



^P5Pe^n 5 n B j 



'i;qpiP2 i]qpiP2 ' 



T 

r,QP3P4 



(48) 



n 5 = n i 1 + E n 2 2 mod 2 

Ql Oil 

n 6 = E "a" + E m ° d 2 ' ( 49 ) 



FIG. 6: Split one rank-four tensor into two rank-three tensors 
and a rank- two tensor. 



We see that n§ and n 6 are completely fixed by n" 1 , / = 
1,2,3,4. Eqn. (48) define three Grassmann tensors (see 
(44)) 



E ^ 



J "j;P6P3P4 / y ' r,PSP3P4 

{™ 6 }{™° 3 }{< 4 } 

^ij-P&Pe ~ (1 T TlSV^fipspn 



(d^rnwT^n^ 4 )" 4 



(50) 



This way, we split the tensor T into to the above three 
tensors, since T can be expressed in terms of the three 
tensors as in (43). It is interesting to note that Gij re- 
duces to a very simple form after one step of second move. 

In the third move, we combine three rank-three tensors 
Ti, Tj, Tk and three rank- two tensors Gij, Gjk, Gu 
into one rank-three tensor T (see Fig. 7). If we 
expand the tensors Ti- Pl p 4 p 9 , j\p2p^p^~! -^^;P3P8P7 &ud 



Gij- P4P5 , Gjk-pepj, G ki .. P8pql we obtain the follow- 
mg coefficients 7"^^ 9 \ T^Jp^ 

T {«3 3 }{«8 8 }{™? 7 } a {«4 4 }{«? 5 } C {«6 6 }{«? 7 } 

k;P3PsP7 a ' y ij:,P4P5 ' y jk\p 6 p 7 

£7fc"p 8 pg" 9 The coefficients of the resulting ten- 

sor' Tj, lP2j , 3 are given by 



77i£> { " 22}{ - 3} = E E E E E E E r ( ^ 8)G -< 9) x 

P4P5P6P7P8P9 {»;*} {n° 5 } {n° 6 } {n° 7 } {n° 8 } {n° 9 } 

r {<H{< 4 }{< 9 }^< 2 }{< 6 }{<n r {™3n{<nK 7 } c {< 4 }{<n a {»6 6 }{»? 7 } c {< 8 }{< 9 } 

i;piP4Psi j;P2P6P5 k;P3P8P7 ^ij;P4P5 ^jk;peP7 =^ki;pspa v ' 
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- T k 



G 




r jk 



FIG. 7: Combine three rank-three tensors and three rank-two 
tensors into one rank-three tensor. 




FIG. 8: A honeycomb lattice and its two sublattices. The 
distance between two nearest neighboring sites is chosen to 
be 1. 



Just as in the bosonic TERG case 24 , the above coarse 
graining transformation can also be generalized onto 
other plainer graphs, such as square, kagome, triangular 
lattice, etc. However, on generic graphs, especially those 
graphs with long range connections, the tensor contrac- 
tion still can be exponentially hard in most cases. 



VIII. EXAMPLE: FREE FERMION STATE ON 
A HONEYCOMB LATTICE 

To test the coarse graining procedure for the Grass- 
mann tensor network discussed above, let us study a sim- 
ple example of a free fermion state and its tensor network 
representations. We assume fermions live on the vertices 
of a honeycomb lattice. Let us consider the pairing state 



l*/> 



'10) 



(52) 



on a honeycomb lattice. Here the pairing operator is 
given by 



E 

(ij) 



uc\c] = 



E < 

ieA,a=l,2,3 



where i and j label sites and (ij) labels the nearest- 
neighbor links of the honeycomb lattice. We know that 
the sites of the honeycomb lattice can be divided into 
two sub-lattice: A and B (see Fig. 8). The pairing of the 
fermions is only between the two different sublattices. 
Here we have used the convention that in the link label 
(ij), i be long to the A-sublattice and j belong to the 
B-sublattice. The three vectors S a , a = 1,2,3, are the 
three vectors that connect a A-site to to its three nearest 
neighboring B-sites. 
Introduce 



c A (k) = 



N 
N 



i£B 



where N is the total number of lattice sites, X^ieA sums 
over the sites in the A-sublattice and J2ieB over tne B- 



sublattice. We can rewrite P as 

p = E c A(- fc H c l( fe ) 



with 



Therefore 



|*/)=n[ 1 + «fe c A(-fe)4(fe)]|0). (53) 



Let us rewrite 



[l + a fe c4(-fc)4(fc)]|0> 
oc [v k c\(-k) + u k c B (k)][u k c A (-k) - v k c f B (k)}\0) 



where 



We find 



\u k \ 2 + K| 2 = 1, u k = v k /u k . 



v k = 



v / T+kF : 



u k = 



v / TTKF 



(54) 



(55) 



(56) 



Let 



ipi-k = [u k c A (-k) - v k c' B (k)], 
-02,fc = [v k c A (-k) +u k c B (k)], 



(57) 



which satisfy the standard commutation relation for 
fermion operators. By construction, we see that 

V>i,fe|%> =V>2,fc|$/> =0. (58) 
Thus |*/) is the ground state of the following quadratic 
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Hamiltonian 

H = + W^l^.k + $l k ifa, k ) 

k 

= ^2l c A(- k ) - a* k c B {k)][c A {-k) - a fe cjg(fc)] 

k 

+ 2KcA(-fe) + cl(k)}[a k c\(-k) + cs(fe)] 
= 2 [ - 2c] 4 (-fc)a fc c| 3 (fc) + h.c. 

k 

+ ^(l ~ \a k \ 2 )[c A (k)c A (k) + 4(fc)c B (fc)] 



(59) 



In real space, the above Hamiltonian can be rewritten as 
H = 2 2|a fe | 2 - 2( 2uc H + /l - c ") 

+ 2(i-3| u | 2 )4 Ci - 2 \ u \ 2 4+^ (60) 

i i,J=l,...,6 

where {A/} are six vectors <5i — ^2, — ^1 + ^2, ^2 — ^3, 
—S2 + ^3, S3 — Si, and —S 3 + Si. Note that, if we 
do a particle-hole conjugation on the B-sublattice, the 
above pairing Hamiltonian becomes the following hop- 
ping Hamiltonian 



H 



'■Cj + h.c. 



(iJ) 

2(1 -3\u\ 2 )(-) i c\c i + Const 



2 R>lM+A^ 

(61) 



i,/=l,...,6 



where (— ) l = 1 if i is in A-sublattice and (— ) l = — 1 if i 
is in B-sublattice. 

Fig. 9 shows the expectation values of the nearest 
paring term c|ct, next nearest neighbor hopping term 

cj +Aj Cj as well as onsite fermion number term c\ci cal- 
culated from the Grassmann-number tensor renormaliza- 
tion algorithm. We found very good agreement with the 
exact results are found for a large range of parameter u. 
The SVD truncation dimension is set to be D cut = 32 
throughout the whole calculation. Fig. 10 shows the 
D C ut dependence for the physical quantities we calcu- 
lated. We find excellent agreement with exact results for 
large enough D cut . 



IX. SUMMARY 

The projective construction has played an very impor- 
tant role in understanding strongly correlated systems. It 
can explain the emergence of fcrmions, anyons (Abelian 
and non- Abelian) , and gauge theory in quantum spin liq- 
uids and quantum Hall states. Mathematically, the pro- 
jective construction also provides an efficient encoding 
for many-body bosonic/fermionic states. 



0.4 



</> 0.3 




FIG. 9: The expectation value of the nearest neighbor par- 
ing, next nearest neighbor hopping and onsite fermion num- 
ber terms from the Grassmann-number tensor-entanglement 
renormalization algorithm. We compare our calculation with 
the exact results and find a good agreement for a large range 
of u. 




30 40 
Dcut 



FIG. 10: The expectation value of the nearest neighbor par- 
ing and onsite fermion number at different D cut (the trunca- 
tion dimension in the SVD decomposition) in the Grassmann- 
number tensor-entanglement renormalization algorithm. It 
shows all those values converge to the exact one at sufficient 
large D cut . In this calculation we choose u = 0.5. 



Recently, Fermionic projected entangled pair states 
have also been introduced as an alternative method 
for efficiently encoding many-body fermionic states. 
In this paper, we show that the strongly correlated 
bosonic/fermionic states obtained from the projective 
construction and fPEPS approach can be represented 
systematically as Grassmann tensor product states. 
These Grassmann tensor product states allow us to en- 
code many-body bosonic/fermionic states with a num- 
ber of parameters that only scales polynomially in the 
number of fermions. We have also shown that it is pos- 
sible to generalize the tensor-entanglement renormaliza- 
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tion group (TERG) method for complex tensor networks 
to Grassmann tensor networks. This allows us to perform 
an approximate calculation of the norm and average local 
operators of Grassmann tensor product states in polyno- 
mial time. 

In conclusion, Grassmann tensor product states can 
be a starting point for a new variational approach to 
strongly correlated bosonic/fermionic systems, as they 
not only include tensor product state and fermion ELF 
states, but also systematically generalize the slave- 
particle projective construction. All the physical prop- 
erties of these states can be efficiently calculated based 
on the Grassmann tensor-entanglement rcnormalization 
group (GTERG) algorithm. Many non-trivial fermionic 
models and frustrated spin models, such as t-J model, 
Hubbard model, as well as Kagome Heisenberg model 
will be studied in our future work. 
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51164 ,the ERC grant QUERG and the FWF grant Fo- 
QuS. XGW is supported by NSF Grant No. DMR- 



0706078. 



Appendix A: Calculate the physical measurements 

In this section, we explain how to use coarse graining 
transformation to calculate the physical quantities of a 
Grassmann tensor product state. As a simple example, 
we first explain how to calculate the expectation value for 
nearest neighbor electron pairing term (\I'|cjcj|\E'). Such 
a term can be represented as a Grassmann tensor network 
with two impurity tensors, see Eq. (41). 

We can apply the coarse graining transformation for 
the uniform part in the same way as we calculate the 
norm. The extra thing we need to know is how to apply 
the coarse graining transformation for the two impurity 
tensors. It turns out we only need to modify the first step, 
where we combine the two rank three impurity tensors 
into a rank four impurity tensor 



E 

P5P6 



i;P5PiP2 j;p6P3P4. "*j;pbP6) 



(Al) 



with 



J 



T' 



i\P5PiP2 



rpll 

j;peP3P4 



" P1P2P3P4. 



E rfZllf^Ei^r nw 1 ) B?1 II« 2 r l 

n? 1 }!™? 2 ) "5 Q i "2 



e T"]i2if }{<4} n(dC) n6Q6 ri( d % 3 )^ 3 n( d ^ 4 ) <4 

{< 6 }{«3 3 }{< 4 } aS 



E 



T 



-///{™? 2 }{«2 2 }{ 
P1P2P3P4 



(A2) 



Because T' and T" contain an odd number of Grassmann numbers in this case, the relationship of the coefficients 
should be modified as: 

r ///{ri2 2 }{< 2 }{n3 3 }{< 4 } = \p \p ( \(J2 a5 n^) T dnt 5 }{n" 1 }{n^} G {n^}{n^} T ,dn^}{n^ 3 }{n^} , 

P1P2P3P4 / j / t \ ' »;P5PlP2 y ij;P5P6 j;P6P3P4 \ ' 

P5P6{„°5 }{„"«} 



The second step is the same as we calculate the norm 
because the rank four impurity tensor contain an even 
number of Grassmann numbers. Of course the last step 
also remains the same except we need to combine one 
new impurity tensor with two new uniform tensors to 
produce a coarse grained impurity tensor, as we do in 
the standard TERG algorithm, see in Fig. 11. If we 
calculate the expectation value for two nearest neighbor 
bosonic operators, such as (niUj), even the first step dose 



not need to be modified because the impurity tensors 
contain even number of Grassmann numbers in this case. 

For generic interactions as well as correlation func- 
tions, the coarse graining procedure can be designed in 
the same way as in the TERG algorithm. Fig. 12 shows 
how to do the coarse graining transformation for generic 
six body interactions on the hexagon of a honeycomb 
lattice. However, we need to take care of the sign factor 
when we apply the coarse graining transformation for 



13 



Y Y 

A A 



i=> 



(a) 



(b) 



FIG. 11: A schematic plot of how to apply coarse graining 
transformations for two nearest neighbor impurity tensors on 
a honeycomb lattice. The filled dots represent impurity ten- 
sors and open dots represent uniform tensors. In (a), we first 
combine two rank three impurity tensors into a rank four 
impurity tensor and then split it into two new rank three im- 
purity tensors. In (b), we apply the last step to combine one 
impurity tensor with two uniform tensors to produce a new 
impurity tensor. 




FIG. 12: A schematic plot of the coarse graining transforma- 
tions for generic six body interactions on the hexagon of a 
honeycomb lattice. The filled dots represent impurity tensors 
and open dots represent uniform tensors. 



impurity tensors which contain an odd number of Grass- 
mann number and in the most generic case, all the three 
steps need to be modified. For example, in the second 
step, if the rank four impurity tensor T'" contain an odd 
number of Grassmann numbers, then we need to decom- 
pose it into two new rank three impurity tensors where 
one has an odd number of Grassmann numbers and the 
other has an even number of Grassmann numbers. We 
also need to take care of the sign factor in the last step if 
we combine impurity tensors with odd number of Grass- 
mann numbers. 



Appendix B: A review of projective construction for 
fractional quantum Hall state 

In this section, we are going to review two examples 
of fractional quantum Hall state to demonstrate how 
to obtain low energy effective theory from a projective 
approach. 8,10 

To use the projective construction to study the v = 1/3 
Laughlin state of N electrons, we consider a system of 
three kinds of partons described by the fermion operators 
ij)a{z), a — 1, 2, 3. The parton system contain N particles 
for each kind of partons, and each kind of parton form 
v = 1 state (which is a ELF state). Let us denote the 
v = 1 state for the a th partons as \^i) a - Then the total 
ground state of the system of three kinds of partons is 
given by ® |Wi) 2 ® ^1)3 = Note that 

the state is a ELF state. Then the v = 1/3 

Laughlin wave function can be expressed as a projection 
of the ELF state 

*3({«i}) = (o| ]l[M*)Mx)M*i)]\*iMi)- ( B1 ) 

i 

Note that the projection simply combine the three par- 
tons into a single electron: 



c(z) = ^1(3)^2(2)^3(2)- 



(B2) 



Although the expression (Bl) is very formal, the pro- 
jection that it describes can be done at the field theory 
level, which allows us to calculate the low energy effective 
theory of the Laughlin state. 

To do the projection at the field theory level, we start 
with the Lagrangian that describes three independent 
partons (that is before the projection) 



Co 



id t ip a 



1 

2m 



igA)Va| 



(B3) 



where we have assumed that each kind of parton carries 
e/3 electric charge and A is the vector potential that 
describes the uniform magnetic field. The chemical po- 
tential is chosen such that each kind of parton fills its first 
Landau level and forms the v = 1 quantum Hall state. 

We note that the above Lagrangian that describes the 
independent partons before the projection has an SU(3) 
symmetry: 



(B4) 



The theory of independent partons contains fluctuations 
of the density and the current of the SU(3) charge. On 
the other hand, we see that the electron operator c trans- 
forms as 



(B5) 



since the i\) a anticommute with each other. Thus the 
electron operator is invariant under the 577(3) transfor- 
mation. As a result, the electronic state obtained after 
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the projection in (Bl) is also invariant under the SU(3) 
transformation since both |0) and c are SU(3) invariant. 
This means that, after the projection, the electronic state 
contains no fluctuations of the density and the current of 
the SU(3) charges. 

This motivates us to perform the projection at the field 
theory level by including a SU(3) gauge theory in the 
above indcpcndent-parton model: 



£ P = 'y]ipti[8abdt - i(ao)ab]-06 



(B6) 



a, b 



a, b 



where (ao,a x , a y ) are the SU(3) gauge fields which are 
3 by 3 hermitian matrix valued fields. The SU (3) gauge 
fields remove all the SU{3) density and current fluctua- 
tions. In other word, if we perform the path integral of 
the SU (3) gauge fields first 



-i / dtd 2 x £„ 



D[a ]e 



- i f dtd 2 x C p 



(B7) 



the resulting effective theory C e will contain no SU(3) 
fluctuations. Thus we say that the path integral of the 
SU(3) gauge fields, / D[a]D[a n } 7 performs the projection 
at the field theory level. 

The full theory is described by the path integral over 
both parton fields ip a and SU (3) gauge fields (ao, a x , a y ): 



/*„]/ 



D[a]D[a ]e 



- i f dtd 2 x C p 



(B8) 



If we exchange the integration order: 



Z = J D[a}D[a ] J D[^ a ] e" 1 **** c » 

= J D[a]D[a }c- i f dtd2xC % (B9) 



we obtain an effective theory that contain only the SU(3) 
gauge field £ a (ao,a x ,a y ). Since each kind of parton 
forms the v = 1 quantum Hall state, the resulting 
SU(3) effective theory turns out to be the level- 1 5?7(3) 
Chern-Simons theory. We conclude that the v = 1/3 
Laughlin state is described by the level-1 5/7(3) Chern- 
Simons topological field theory. All the topological prop- 
erties of the v = 1/3 Laughlin state, such as fractional 
charges and fractional statistics can be obtained from 
such a SU (3)i Chern-Simons theory. Note that the level- 
1 SU(3) Chern-Simons theory is equivalent to a U(l) 
Chern-Simons theory. Thus SU(3)\ Chern-Simons the- 
ory actually describes an Abelian state. 

We have seen that if we let partons to form the v = 1 
quantum Hall state, we will obtain the v = 1/3 Laughlin 
state and the SU(3)i Chern-Simons theory as its low 
energy effective theory. If we let each of the three kinds 
of partons forms the v = to quantum Hall state, we will 
obtain a v = m/3 quantum Hall state and the level- to 



SU(3) Chern-Simons theory as its low energy effective 
theory. 10 Such a v — m/3 quantum Hall state is a non- 
Abelian quantum Hall state. 

Similarly, starting from the ELF state for four kinds of 
partons | 1 3/i\E f i\E'i 1 3>i), we can construct the non- Abelian 
v = 1 Pfaffian state (for bosonic electrons) 

*p&({*}) (BIO) 

= <o| Y[[M*i)M*i) - ^3(«i)ife(*i)]|*i*i*i*i>- 

i 

where the electron operator is related to the parton op- 
erators as 



c(z) = xp!(z)xp 4 (z) - ip 3 (z)ip 2 {z). 



(Bll) 



To do the projection at the field theory level, we start 
with the independent parton model 



4 

C = [^i id f<l>a 



a=l 



(B12) 



We note that the above Lagrangian has a SU (4) symme- 
try: 



V'a ->• U ab ^ U G SU(4). 



(B13) 



Thus the theory of independent partons contains SU(4) 
charge and current fluctuations. 

However, the electron operator c = tpi^i ~ ^3 ^2 is not 
invariant under the full SU(4) transformations. If we 
identify (ipi, ipt) = (ipu, V>i2, ip2i, ^22), we find 

c = ^11^22 - ^21^12 ( B14 ) 
oc ^11^22 - ^21^12 + ^12^21 - i>22i>n 

OC lp aa (T2)ab(<?l) a ^bl3 OC 1p T \t 2 ® U 2 ){t <g> C7 3 )^, 

where r = cr are the 2 by 2 identity matrix and r/, (J; 
are the Pauli matrices. Here 77 acts on the first subscript 
a of ip aa while 07 acts on the second subscript a of ip aa . 
We see that the electron operator is invariant under a 
subgroup of SU(4) generated by 10 generators: n ® ao, 
t,®<7i, Ti®T2, and tq®(J3, It turns out that the above 10 
generators generate the SO(5) group in its 4 dimensional 
spinor representation. Therefore, the electronic states do 
not contain any SO(5) fluctuations. 

To remove the 50(5) fluctuations at the field theory 
level, we can include a SO (5) gauge field (the 4 dimen- 
sional representation spanned by Tj ® ao, r» ® a\, Ti®U2, 
and ro <8> 03) in the parton Lagrangian. After integrating 
out the partons we obtain the low energy effective Chern- 
Simons theory for the Pfaffian state, which is a SO(5) 
Chern-Simons theory. A different effective Chern-Simons 
theory for the PfafRan state was obtained in Ref. 28. 

We would like to point out that in the first form of the 
projective construction (Bl), the electron operator c(z) 
is expressed as a product of parton operators. In this 
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case, we can use variational Monte Carlo calculation to 
numerically study many properties (such as ground state 
energy) of such projective states. On the other hand, in 
the second form of the projective construction (BIO), the 
electron operator c(z) is expressed as a sum of several 
products of parton operators. In this case, in general, 
the variational Monte Carlo method is ineffective due to 
the sign problem. So far, we still do not have an effective 
numerical method for the second form of the projective 



construction. 

The Grassmann tensor network can represent projec- 
tive states obtained from the both forms of the projective 
construction. So the renormalization of the Grassmann 
tensor network might allow us to approximately calculate 
the norms and average local operators for both forms of 
projective states. A detailed study will be presented in 
our future work. 
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